Phonon-assisted tunneling in the quantum regime of Mni2-ac 
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Longitudinal or transverse magnetic fields applied on a crystal of Mni2-ac allows to observe inde- 
pendent tunnel transitions between m=-S+p and m=S-n-p (n=6-10, p=0-2 in longitudinal field and 
n=p=0 in transverse field). We observe a smooth transition (in longitudinal) from coherent ground- 
state to thermally activated tunneling. Furthermore two ground-state relaxation regimes showing 
a crossover between quantum spin relaxation far from equilibrium and near equilibrium, when the 
environment destroys multimolecule correlations. Finally, we stress that the complete Hamiltonian 
of Mni2 should contain odd spin operators of low order. 



75.50.Xx, 75.45.+j, 71.70.-d 

The system Mni2-ac is constituted of magnetic 
molecules in each of which 12 Mn ions (8 spins-2 and 
4 spins-3/2) are strongly coupled by super-exchange in- 
teractions pL The resulting spin S — 10, containing 
10 3 atoms of different species (Mn, O, C, H) has volume 
sa lnm 3 . The Hilbert space dimension Djj describing the 
entangled spin states of the molecule is huge: Djj = 10 8 . 
In the mesoscopic approach to Mni2-ac most observa- 
tions can be explained with Dh = 2S + 1 = 21 (for a 
review, see [||). The large energy barrier DS^ gives ex- 
tremely small ground-state splitting. The consequence 
is that quantum relaxation is slow and strongly influ- 
enced by fluctuations of the environment. Above 1.5 K 
a single crystal of Mni2-ac shows staircase- like hysteresis 
loops when the magnetic field is applied along the c— easy 
axis of magnetization j^] (for powder experiments see 
[0-tfJ )■ Tunneling rates, deduced from relaxation and ac- 
susceptibility experiments give sharp maxima precisely at 
the fields where hysteresis loops are steep: B n « 0.44n 
(in Teslas), with n = 0, 1, 2, . . . All these results consti- 
tuted the first evidence for thermally assisted tunneling 
in large spin molecules. Later we studied the effects of 
spin and phonon baths in this tunneling regime and the 
passage to the classically activated regime ( [|]|| and |Tq1 
for a low spin molecule). 

In this letter we give a similar study, but at lower 
temperatures. We determine the nature of the cross- 
over between thermally assisted and ground-state tun- 
neling and show an interesting quantum behaviour of the 
molecule spin dynamics, related to spin-phonon transi- 
tions. Furthcrcmore, we suggest that odd spin opera- 
tors pH , consistent with the tetragonal S4 point sym- 
metry of Mni2-ac |j| should be taken into account in 
the Hamiltonian of this system. The magnetization of a 
small single crystal of Mni2-ac in high fields and at sub- 
Kelvin temperatures was obtained from magnetic torque 
experiments, performed at the Grenoble High Magnetic 
Field Laboratory. The way to extract the longitudinal 
or transverse magnetization (parallel or perpendicular to 
the easy c— axis) from the measured torque T is analyt- 



ical and unambiguous. The applied magnetic field cre- 
ates a torque, which at equilibrium, is compensated by 
a mechanical drawback one. In a longitudinal field the 
system of non-linear equations linearizes. To find a hys- 
teresis loop we first perform a two-parameter linear fit 
(u and v) of the saturation region with Bq/T = u + vBq 
and then we calculate the normalized magnetic moment 
to = M\\/Mg — u/{Bq/T — vBq). In the transverse 
case, M± is, to a good enough approximation, given by 
M X = T/B . 

Fig. [l] shows several staircase- like hysteresis loops in a 
longitudinal field, below 1.3 K. The angle between the 
applied field and the easy axis, due to some inherent 
misalignment in the experimental set-up, was < 2.5°. 
The corresponding transverse field component, smaller 
than 0.2 T, is almost not affecting the transitions with 
2to — n a multiple of 4 (to is the S z eigenvalue) , because 
of the leading C term in the transverse anisotropy (see 
the caption of Fig. |]b) or with n < 6. The observed B n 
are independent of temperature (Fig. |^), an effect which 
was also observed in but below « 0.7 K only. This 
suggests that tunneling takes place from the metastable 
state m = —10 to final states to = 10 — n (n = 6 to 
10). As shown in Fig. [j] inset for n = 8, the B n are de- 
termined accurately enough to see that each resonance 
is constituted of several temperature-independent peaks 
(resonance groups). Their widths « 0.15 T are close to 
the one observed at high temperature and come from the 
distribution of dipolar fields . A comparison between 
the measured and calculated B n makes their interpreta- 
tion obvious. The calculated B n are given by level cross- 
ings in Fig. ||b (the Hamiltonian, given in the caption 
of Fig. |b, was found by EPR experiments jL3)). Inter- 
cepts of the to = S = — 10 state with m = 3, 2, 1 and 0, 
give the B n indicated by horizontal lines n-0 in Fig. ||a. 
Similar intercepts of the first excited level to = —9 with 
to = 2, 1, and -1, are also given. This shows unambigu- 
ously that tunneling occurs, when increasing tempera- 
ture, first from m = —10, then from the first excited state 
to = —9, then from m = —8, . . . Tunneling transitions 
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with smaller m are closer and closer and may appear to be 
temperature-dependent if the field resolution is not good 
enough. Regarding peak intensities, the area of decon- 
volved Lorentzian peaks in each resonance group, char- 
acterized by the same index n, determines the fraction 
of tunneling events taking place from the ground-state 
(index p = 0), from the first excited state (p = 1), ... 
One can define an experimental cross-over temperature 
T Cj „ for which the two peaks (n-0 and n-1) have the same 
area; for example, for n — 8 we found 0.95 K. Account- 
ing for the proportionality of the tunneling rate with A 2 , 
where A is the tunnel splitting |H , this crossover can be 
written kT c , n w S ,i, n /21n(A n _i/A n _o), where £ ,i, n is 
the separation between the ground-state and the first ex- 
cited state, near the n th crossing. For the resonances 
with n = 8 and £?o,i,8 ~ 7.3 K, calculated by exact diag- 
onalization, we get ln(A 8 _i/A 8 _ ) = 3.88 corresponding 
to a transverse field of « 0.31 T (a misalignement of 
« 4.5°). Obviously, such a transverse field is larger than 
what is expected from our experimental set-up. Simi- 
lar discrepencies are also observed with other resonances. 
The important effect of small transverse fields is not suf- 
ficient to give an overall agreement, in particular when 
2m — n is not a multiple of 4 and when n > 6. Finally, 
early observation of a crossover temperature of w 2 K 
in zero field P,p|, p^[ is now confirmed by the extrapo- 
lation of T C) „ (Fig. |2|a), giving T Ci0 ~ 1.7 K (the same 
type of extrapolation shows that the continuous process 
of ground-state to activated tunneling up to the top of 
the barrier spreads between w 1.7 and «s 2.4 K). How- 
ever, this requires a ground-state splitting w 3 x 10~ 10 K 
instead of the w 10~ n K as calculated in This zero- 
field crossover All that suggests a lack of knowledge of 
the transverse spin operator coefficients of the Hamilto- 
nian used for Mni2-ac (Fig. 2b). Besides the fact that 
all crystal field parameters of order larger than four are 
unknown, we must mention that all operators with odd 
order are missing. The S4 point symmetry of Mni2-ac 0] 
leads to the following Hamiltonian 0: H = B^O% + 
BlO\ + B\0\ + B\0\ + Bl0 2 5 + B° 6 0° 6 + B|0| + . . . where 
the 0™ are the Steven's equivalent operators Jl^| . The 
complete Hamiltonian should also contain magnetic field 
and the effect of Dzyaloshinsky-Moriya terms |l5|,[l7],|l8| . 
Interestingly, the new odd terms have no diagonal con- 
tribution, which is consistent with the good agreement 
we obtain regarding the resonance positions. This also 
shows that diagonal terms of order > 6 are negligible. 

Let us now consider the case when the field is applied 
perpendicularly to the easy axis. The symmetry of the 
double well is preserved allowing to observe tunneling 
between the ground-states \S, +S > and \S, — S >. The 
sample was first saturated along the c— axis of the crystal. 
Then, the field was set to zero and the sample rotated to 
90° ±1°. During the experiments, the applied field Bq 
was varied between and 5 T with a transverse compo- 
nent Bt nearly equal to Bq, and a longitudinal one Bt 



varying between and « -0.07 T. When B L w -0.04 T 
w the Lorentz internal field of Mni2-ac, the internal field 
is canceled. Nevertheless, level widths being « 0.15 T, 
the spin-up and spin-down density of states were in co- 
incidence during the whole variation of the applied field. 
Hysteresis loops of the transverse magnetization are plot- 
ted Fig.^. They show only one jump and are nearly 
independent of temperature, indicating that spin rever- 
sals take place by single coherent ground-state tunneling 
events. This loop can be understood easily. In trans- 
verse fields smaller than w 4 T, tunneling rate Tq between 
symmetrical ground-states is too small and nothing hap- 
pens. The tunneling gap increasing almost exponentially 
with Bt jn|) T becomes rapidly very large, and when 
Bt ~ 4.5 T the spins tunnel to the other well during the 
time of the experiment; due to a small misalignment from 
the ideal transverse geometry, as shown in Fig. |^ insert, 
after tunneling the molecules are kept within the final 
well (the sample is again saturated) . The observation of 
tunneling requires T^ 1 to fall in the experimental win- 
dow 1 — 10 3 s. Such values for the tunneling rate, which 
depends on the level splitting and the hyperfine field fluc- 
tuations H], are realistic according to NMR results on 
Mni2 Q and for tunnel splittings calculated between 4 
and 4.5 T. 

Magnetic relaxation experiments give direct access to 
tunneling rates. An example is given Fig. ^, in longitudi- 
nal and transverse applied fields. In both cases and be- 
low 1.2 K, the relaxation is a square root of time at short 
times: M(0) - M(t) = M sy /T sqrt t. This law, already 
observed in Fes |2l[] or Mni2-ac above 1.5 K S, was pre- 
dicted by Prokof 'ev and Stamp jl4j . It is characteristic of 
quantum tunneling in the presence of dipolar interactions 
and fast hyperfine fluctuations and ends for times « rj grt 
due to the building of multimolecule spin correlations. 
At longer times more and more important correlations, 
should in principle lead to non-exponential relaxation 
(see Fig. 1 in fl4| ). However, this is not what we observe 
in Fig. where the quantum relaxation regime, succeed- 
ing the square root one at long times, is exponential. 
We believe that this is because during the smooth tran- 
sition between the square root and exponential regimes, 
couplings to the environment (spins, phonons. . .) have 
enough time to erase multimolecule spin correlations, an 
effect which could be called "environmental annealing". 
The crossover time t c (when the dashed fit intercepts 
the dotted fit in Fig. ^), depends on many parameters 
and will be analysed in further studies. Regarding the 
first regime and for our present experimental conditions 
(sample in good contact with the cryostat), Tr^T) -1 
is always of the order of t c , showing that fast recovery 
has the same origin as fast spin relaxation, i.e. not too 
small tunnel splittings and large hyperfine field fluctu- 
ations Jl4| , p2[ . As shown above spin-phonon transitions 
also contribute to Tt,l and thus to recovery near T Cj „ 
by adding new, thermally excitated channels, with larger 
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tunnel splittings and faster relaxation rates. 

In Fig. | inset are shown the square-root relaxation 
rates plotted vs. temperature, measured at the 8-0 reso- 
nance (Fig. ||) in longitudinal geometry (full circles-r^) 
and in transverse geometry for the 0-0 resonance (full 
squares-Fr). As expected from the first part of this pa- 
per, Tl increases above w 0.6 K, then it goes through 
a maximum and then decreases. The increase of Tl(T) 
is due to the proximity of the resonance 8-1. The field 
separation between this resonance at 3.90 T and the res- 
onance 8-0 at which Tl(T) is measured at 4.02 T, being 
of the order of level broadening by dipolar interactions, 
spin-phonon transitions between them are allowed. This 
is what we observe and constitutes a clear example of 
phonon-assisted tunneling in the quantum regime. Re- 
garding peak intensities, near equilibrium the rate in- 
creases like T L (T) « r L:8 _ + r i)8 -i exp(-E ,i, s /kT). 
As an example, for T = 0.7 K and Fl i8 _ = Ti(0.55), we 
obtain ln(A 8 _i/A 8 _o) ~ 5.35 too large to be explained 
numerically. As we already mentioned, new terms in the 
spin Hamiltonian of Mni2 could lead to a better agree- 
ment. The decreasing part of Tl(T) simply comes from a 
decreasing of the number of states available for phonon- 
assisted tunneling, due to the magnetization decrease 
when the field sweeps across the 8-1 resonance. Like Tl, 
the square-root relaxation rate r^, measured in trans- 
verse geometry, is independent of temperature at low- 
est temperatures, as expected for ground state tunneling. 
Above w 0.7 K the tunneling rate increases with temper- 
ature. However this is a very smooth effect which cannot 
be due to activated tunneling on levels to — 9 and to = 8. 
It is rather due to tunneling assisted by spin-phonon tran- 
sitions between m w ±10 states, split by dipolar fields 
(and integrated over the spin-up and spin-down distribu- 
tions) . This assumption is supported by the fact that the 
corresponding energy scale of ~ 0.5 K for increasing Tt 
is close to the dipolar energy (g^sS 'B dipolar ~ 0.5 K). 
This is another example of phonon-assisted tunneling in 
the quantum regime. 

In conclusion, we have shown that in Mni2-ac with a 
longitudinal field (tj ^ 0), the transition between coher- 
ent tunneling on the ground-state and on upper states 
(with thermal activation) takes place, level after level, 
when they come into coincidence in an applied sweep- 
ing field. This is an obvious consequence of the 4 t/l or- 
der anisotropy term BS^, which prevents the spin-up 
and spin-down level schemes coinciding. In any case, 
the experiments show a continuous transition, even in a 
constant field, contrary to Ref. |^|. More interestingly, 
phonon-assisted tunneling is evidenced between m = — 10 
and to — 10 splitted by dipolar fieds in the transverse 
case or between m = — 10 and to = 1 in the longitudi- 
nal one. "Annealing" of multimolecule spin correlations 
in the quantum regime result from spin relaxation itself 
and from spin-phonon transitions near the crossover to 
excited states. Such an exponential phonon-mediated re- 



laxation cross-over was predicted for a single-molecule re- 
laxation p^ ], but no theoretical work exists which takes 
into account dipolar interactions. It would be of great 
interest to see a theoretical treatement of the long time 
quantum relaxation in the presence of phonons and dipo- 
lar interactions. Finally we suggest not to forget that 
spin operators with odd powers are allowed by the S4 
symmetry of the Mni2-ac molecule. 

We thank I.Tupitsyn and S. Miyashita, S. Maegawa 
and T. Goto for discussions. 
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FIG. 1. Hysteresis loops of Mni2-ac, derived from torque 
measurements with an applied field along the c— axis. Below 
ps 0.7 K, the curves are temperature-independent, indicating 
quantum tunneling without thermal activation. Transitions 
widths result from the distribution of dipolar fields. Inset: 
dm/dBo where m — M«/M 3 for the resonance n — 8 at three 
temperatures: 0.90 K (dots), 0.95 K (dashed) and 1.00 K 
(continuous). The resonance is split in two: tunneling from 
the ground-state m = —10 and from the first excited state 
m = 9. Note the change of the height of each peak, when 
varying the temperature. Similar behavior was found for all 
the other resonances. 




FIG. 2. (a): Measured positions of peak maxima 
vs. temperature. Horizontal lines indicate the cal- 
culated crossing fields. (b): Energy levels spectrum 
calculated by exact diagonalization of the Hamiltonian 
H = -DSl - BS A Z - C{S\ + Si) - g^ B S z H z for a spin 
S = 10 and the parameters D — 0.56 K, B = 1.1 mK, 
C = 3 x 10~ 2 mK obtained in EPR. The term B shifts the 
resonances belonging to same n and different m to lower fields 
when m decreases. 
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FIG. 3. Half hysteresis loop in transverse field at three dif- 
ferent temperatures. M± is given by the measured torque 
divided by field. The crystal was first saturated in a large 
negative field. Inset: the geometry of the experiment. 
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FIG. 4. Relaxation curves in a transverse (longitudinal) 
field of 4.42 T (4.02 T) for the resonance n = at T = 0.5 K 
(n = 8 at T = 0.9 K). In both cases the square root regime 
(dashed) is observed only at short timescales, then the relax- 
ation becomes exponential (dotted). Note that experiments 
are done far from saturation. Inset: temperature dependence 
of the square root tunneling rate in transverse (squares-Fr) 
and longitudinal field (circles-Pi). 
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